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Abstract 

This report documents the theory and experimental work which demonstrated 
the feasibility of conical-scan tracking with the NASA/ Jet Propulsion Laboratory 
64-m-diameter paraboloid antenna at Goldstone, California. The purpose of this 
scheme is to actively track spacecraft and radio sources continuously with an 
accuracy superior to that obtained by manual correction of the computer-driven 
pointing. The conical-scan implementation gives increased tracking accuracy with 
X-band spacecraft signals, as demonstrated in the Mariner Venus/Mercury 1973 
mission. ALo, the high accuracy and ease of measurement with conical-scan 
tracking allow evaluation of systematic and random antenna tracking errors. 
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Conical-Scan Tracking With the 64-m-diameter 
Antenna at Goldstone 


I. Introduction 

The 64-m-diameter NASA/JPL paraboloid antenna at 
Goldstone. California, is used for communication with 
many spacecraft. It is also used for research in such 
supporting areas as planetary radar, interferometry, and 
radiometry. For conventional pointing, the antenna 
employs a “master equatorial” mechanical reference 
which is located at the intersection of the azimuth and 
elevation axes. The master equatorial is driven by a com- 
puter, and the antenna is slaved to the master equatorial 
by a large servo system which constitutes the azimuth and 
elevation drive mechanisms. Although high off the 
ground, the master equatorial is atop its own separate 
structural tower and foundation and is protected from the 
wind and sun by appropriate shielding. 

For S-band frequencies in the range of 2100 to 2400 
MH?, the antenna has excellent pointing performance. 
The absolute pointing error is a small fraction of the 
antenna beamwidth at S-band, and for use there the 
antenna “tracks” a spacecraft or radio source by si^.ply 
pointing in a computer-calculated direction. Active track- 
ing schemes using monopulse or conventional conical- 
scan tracking were not installed on the antenna due to its 
excellent absolute pointing ability. 


The Mariner Venus/ Mercury 1973 mission, the first 
mission on which conical-scan tracking was used to track 
a spacecraft, had an X-band downlink as well as the usual 
uplink and downlink at S-band. With approximately a 
factor of 4 increase in frequency over S-band, the antenna 
beamwidth was reduced by the same factor. In mission 
planning, concern arose over the ability to keep the X- 
band antenna beam (approximately 0,038 deg wide at the 
half-power points) properly pointed so a gain loss would 
not occur. Implementation of a monopulse or conven- 
tional conical-scan tracking feed would somewhat degrade 
the system operating temperature and would entail sub- 
stantial hardware modifications. Implementation of a con- 
ical-scan tracking scheme for radio sources by rotating an 
entire 26-m-diameter paraboloid antenna about tracking 
boresight had been done previously (Ref. I). The idea was 
felt to be sound for spacecraft as well, and it was decided 
to implement it on the 64-m-dii:meter antenna. Also, a full 
analysis of the scheme was to he carried out. This report 
documents the results of this effort. 

We will treat angp tracking of a noncoherent radio 
source first in Section III and then, m Section IV, con- 
sider tracking of a coherent source (spacecraft). The 
approach here assumes a basic knowledge of conical-scan 
tracking such as presented by Skolnik (Ref. 2). In Sections 
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III and IV, system gain fluctuations are ignored. This 
effect can dominate the system performance if the con- 
ical-scan rate is slower than the gain fluctuations. The 
reader should be aware that this will probably be the case 
for lobing done by moving the entire antenna. Then the 
results in Sections 111 and IV are optimistic and do not 
give the complete picture. If a faster scanning technique is 
employed in the future (e.g., feed or subreflector nutation) 
we might be able to ignore the gain fluctuations, and then 
Sections III and IV will give the complete syrtem perfor- 
mance. In Section V, the effects of gain fluctuations will 
be considered in detail and the performance for arbitrary 
scan rates will be obtained. To summarize: 

(1) Sections III and IV give tracking performance 
under the assumption that the scan rate is much 
faster than gain fluctuations. 

(2) When lobing the entire antenna, this assumption is 
probably unjustified, but it might be justified if the 
scan rate were significantly increased in the future. 

(3) Gain fluctuation effects are considered in Section 
V. 

The remainder of this report covers the mechanics of 
implementation on the 64-m-diameter antenna (DSS-14) 
and the experimental results obtained thereby. 

II. Preliminaries 

The conical-scan scheme we shall employ here is essen- 
tially that which was proposed by G. Levy (Ref. 3) and 
successfully demonstrated with radio sources by R. Gos- 
line on the NASA/JPL 26-m-diameter paraboloid DSS-13 
at Goldstone, California (Ref. 1). This method had not 
been used with spacecraft until the present investigation. 

Angle tracking is accomplished by lobing the antenna 
around boresight in a circular pattern with constant 
angular offset, called the scan radius or “squint” angle 
(see Fig. 1). “Boresight” refers to the entire system and 
not to the peak of the antenna gain pattern. Pointing 
“predicts” are >ed to direct boresight very close to the 
target and then the conical-scan lobing around boresight 
generates “offsets” to correct boresight for pointing errors. 
Generation of these offsets will be discussed shortly. 

We henceforth assume the antenna beam to be circu- 
larly symmetric because asymmetry of the beam causes 
only second-order effects. If a spacecraft or radio source 
is on boresight, the received signal power in the antenna is 
constant with time. However, if the source is off boresight, 
a small sinusoidal variation in received power occurs. The 
frequency of the sinusoid is simply the conical-scan rate. 



the sinusoid amplitude is proportional to angular devia- 
tion of boresight from the source, and the phase of the 
sinusoid with respect to the conical-scan rotation gives the 
direction of the error. 

Conical-scan schemes for radar tracking use sinusoid 
and cosinusoid signals referenced to the antenna lobing to 
do a correlation on the received signal strength to derive 
error signals for two orthogonal coordinates of the an- 
tenna, say elevation and cross-elevation (i.e., azimuth with 
secant correction). Radar tracking conventionally uses a 
relatively high scan rate of many cycles per second, so 
analog correlation suffices. However, for scan periods of 
10-100 seconds the double-frequency terms arising in 
analog correlation would cause unwanted oscillation in 
the tracking system. Therefore, a scheme using integration 
over exactly one scan period and subsequent offset 
generation was chosen. This scheme is also compatible 
with current pointing offset procedures in the NASA/JPL 
Deep Space Network (DSN). The system block diagram is 
given in Fig. 2. The signal strength in Fig. 2 is taken as the 
output of a suitable radiometer when tracking a radio 
source. Careful attention must be paid to the time con- 
stant of the radiometer as we shall see later. The signal 
strength will be the receiver automatic gain control (AGC) 
voltage when tracking a spacecraft with coherent down- 
link carrier. 

Initially, we shall proceed under the implicit assumption 
that all digital sampling and arithmetic operations are 
perfect; i.e., we ignore the effects of sampling rate and 
quantization noise. We defer consideration of these very 
important problems to a later time. 

We have assumed the antenna beam to be circularly 
symmetric. We define <7(/?) to be the antenna power gain 
at angle ft off peak gain. We do not yet specify the form 
of Gift). In Fig. 3 is depicted the conical-scan geometry. 
We initially assume that the antenna boresight is B radians 
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Rg. 2. Syttoffii block diagram 


away from the source in the cross-elevation axis and <J> 
radians away from the source in the elevation axis. The 
conical-scan pattern is also shown in Fig. 3 as a circle of 
radius R, henceforth called the scan radius. 

We will now calculate the error signal which the angu- 
lar error generates. The instantaneous antenna 

(peak gain) pointing position is given by 

( 1 ) 


4 >' ELEVATION 



Fig. 3. Conical' scan gaomatry 


<(>'a=<j>-h/?sinto^/ (2) 

where is the conical-scan rate in radians/second. This 
scanning is generated by the computer. 

We assume the source to be, at the origin, 9'^0^<f>\ 
Therefore, the instantaneous angle by which the antenna’s 
peak gain is away from the source is, by the Pythagorean 
theorem, for small 9' and </>', 

/3=V(0')^ + («J>')' (3) 

and, by Eqs. (1) nd (2), 

COStJ^/ + <j.sinu„/) (4) 

We assume that 9 and (the errors) are much smaller 
than R since this will be the case in practice and we can 
drop the and terms. 

We thus obtain 

P * }/r^ + 2R9co%o:j + 2/?<>sin (5) 

Now again since 9 and <> are we have R9<t:R^ and 
/?<;»</?, and since |cosu^/| and |sinw^/| are <1, we can 
use 

/nr + i 
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and have 


^ — + + (6) 

Now the instantaneous gain in the direction of the target 
as a function of time G(/) is GIAO) so by expanding 
G(P) in a Taylor’s Series around we have 

G( /) - G( /? ) -f G'( /? )^ cos w J 


ANTENNA 



RADI0^ 

OISPLA'. 


SIGNAL STRENGTH 


G'(/^)<J>sinw^/ 


where 


G'(/?)- 






0 ) 


( 8 ) 


For simplicity we fake G\R) as a positive number when 
R>0 even though it is actually negative. We have thus 
obtained an expression for the system gain in the target 
direction as a furiCtion of time in terms of the angular 
errors 9 and <}> in orthogonal coordinates. Note that when 
that G(r) is constant as one expects when lobing 
around the target at constant offset. 


Fig. 4. Radlomattr configuration 

conical-scan period so that little amplitude and phase 
degradation of the error signal sinusoid results. If the 
radiometer smoothing lowpass filter (LPF) has a time 
constant satisfying this requirement, then the signal 
strength can be taken at the normal radiometer output. 
However, if the radiometer smoothing time approaches 
the conical-scan period, a fast RC filter mu^t be used or 
severe performance degradation will result because the 
sinusoidal error signal will be smoothed out. In Chapter 
IV, the filtering caused by AGC is examined in detail. The 
results there are applicable here as well if the fast RC is 
not fast enough. 

From Eqs. (9) and (7) the system total (inrtantancous) 
temperature input is 


III. Radio Source Tracking 

We now wish to consider the performance of the con- 
ical-scan tracker on a point radio source. We assume a 
conventional total-power radiometer following the an- 
tenna as in Fig, 4. We take the system operating tempera- 
ture to be and take to be the radio source 
temperature when seen with the beam peak gain (i.e., 
conventional radiometry without conical-scanning). The 
radio source temperature seen when the source is on the 
conical-scan borcsighl is clearly 




(10a) 


® ^OF ^ ) 

+ 7’v^X^)[^cosw^/ + </)S Uij] (10b) 

where we have now introduced the normalized gain 


G(/) 

F 's* - r 
G(0) 


( 9 ) 


g(P)^G(P)/G(^) 


and its derivative 


where G(0) is the peak gain of the beam. In Fig. 4 we 
lump all RF and IF filtering into the voltage transfer 
function //(/). We follow this by the conventional 
square-law detector and a fast lowpass resistor-capacitor 
(RC) filter. At this point we take the signal strength to the 
conical-scan circuitry (actually the computer). After the 
fast RC filter we also allow a conventional smoothing 
filter for radiometer smoothing. The fast RC filter is 
chosen to have a time constant no more than I /20 of the 


g\ft)~G\P)/G{Q) 


The output of the square-law detector is 
x(/)»CT,B + w(t) 

where B is the effective bandwidth given by 


(II) 


k 

\ 
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where is power gain and C is just a constant 

representing the receiver gain (Ref. 4). The random pro- 
cess /»(/) has effectively a (two-sided) white power 
spectrum in walls/ Hz given by 

SSf)~C^T.^B^C^[To^ + T^R)YB (13) 

where 0 and arc assumed zero in Eq, (10b) to get this 
last result since they are very small and C is the same as 
in Eq. (1 1). We may assume /i(/) to be a Gaussian random 
proccs:. 

Using Eqs. (10b) and (1 1) we have 

x(0-O[ 

^OP + Tsg(R)] 

+ CBTsg'{R)[ecosu„! + <j>sin w„/ ] + «(/) ( 14) 

To develop the elevation axis offset for each scan we form 


and 

+ (18) 

where we have subscripted 8 and ^ to indicate their values 
during the scan. Also, 


A^CBT^{R)P/2 

(19) 

f^n{t)sinujdl 

•'0 

(20) 

f'’n(i)cosu)jdl 

•'o 

(21) 


The parameter A is readily measured by observing I’he 
result of the integrals in Eqs. (15) and (16) when the 
system (in open-loop operation) has deliberate errors 9 
and <>. 

In the Appendix it is shown that and arc 
zero-mean random variables with common variance 



x(r)sint.:^rd/ (15) 

•'o 

where /i is a selectable gain, the e superscript on A denotes 
elevation axis, and the subscript k denotes the offset that 
is generated by the /c*** scan. We have also defined the 
conical-scan period 

We similarly develop the cross-elevation offset as 

x{t)cosiC„i dl (16) 

-'o 

We choose the same gain h for both axes, as there seems 
to be no advantage in having different gains. Because we 
arc integrating over exactly one period of periodic func- 
tions it is easy to see that 

+ (* 7 ) 


oi- Var[ - Var[ - C^[ To, + T^{R)yBP/2 

( 22 ) 

A. Dynamics of tho Trackor 

We now consider the dynamic performance of the con- 
ical-scan tracker. The borcsight position of the tracker has 
been denoted by and where the index k indicates 
their values during the k^*' scan. We have assumed the 
target to be at the origin. The dynamical equations then 
lake the difference equation form 



(elevation) 

(23) 


(cross-elevation) 

(24) 


where and A^^^ arc the offsets calculated above, which 
arc generated by the k'^ scan. We now will consider the 
elevation d> namics only and will later apply the results to 
the cr^ss-clevation case as well. From Eqs. (17) and (22) 
we sec that the mean and variance uf A^'^ are 
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(26) 


Fiom Eq. (23) and Using Eqs. (25) and (26) wc can 
write 


can be used only when the system SNR is high enough so 
that the error can be accurately corrected in one scan. 

Wc shall hencefc th assume that 0<r< 1 unless speci- 
fied otherwise. When r is near unity (say 0.7-0.9) the 


(27) 


(28) 


where we define the decay factor r as 
r^X-hA 
and 

£[«</’] -0 
Var[«i'*] » Var[A)'>] 


The equation in (27) is a first-order linear difference 
equation with constant coefficients. As such, wc could 
obtain a complete solution but only two characteristics 
arc needed here: (1) response time of the tracker (note 
that it will usually be much longer than the scan period), 
and (2) steady state tracking error. 

Ignoring the solution of Eq. (23) for an initial of^ n 
<f), is 




(31) 


OL 

o 


o 

z 

< 


(a) r = 1 .2 

2r 


2P 3P 

TIM£ t 


4P 


(b) r = 1 

I 


01 

(29) 

U 

0 

(c). 

(30) ,l 

0.7 

r=^ 


2P 3P 

TirAE t 


4P 






2P 


TIME 


iap 4P 

2.81 P 


(d) r - 0 

ir~ 


2P 3P 

TIME t 


4P 




-,l! 


2T 


3P»“ 




5P 


5P 




5P 


5P 


5P 


TIME ♦ 


The decay factor r deternmnes the speed and nature of the 
response of the closed-loop tracker. In Table I is 
summarized the behavior of the conical-scan tracker for 
various values of r. We have assumed that the noise term 
in Eq. (23), is zero, but when the system signal-lo- 
noise ratio (SNR) is high the actual system behavior will 
be very similar to this noiseless case here. In Figs. 5(a) to 
5(g) arc depicted typical responses for each case in Table 
1, where we assume an error of one unit to exist during 
the first scan. From the figures it is clear that cases a, 6, /, 
and g arc of no value since the error does not diminish. 
Although cases c and e both can converge absolutely at 
the same rate if the absolute values of r are equal, we will 
see later that case c will have smaller error due to noise 
than will case e. Therefore, case e is of no value and wc 
arc left only with cases c and d as being of use. Case d 


(f) r -1 







”~1 

0 

p 

2P 

TIME 1 

3P 

4P 

51 


Fig. S. Tracking raaponaaa 
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Tabto 1. Tracking rctponta characltrlttlct 


Case 

Value of r 

Noiseless response 

Example 

a 

r> + 1 

Monotonic divergence 

1-ig. 5a 

h 


Constant value 

Fig. 5b 

c 

0< r< + I 

Monotonic convergence 
(useful case) 

Fig. 5c 

d 

r-0 

Converges in one scan 
(sometimes useful) 

Fig. 5d 

e 

- 1 <r<0 

Oscillatory convergence 

Fig.5e 

/ 


Constant amplitude oscillation 

Fig. 5f 


r<-\ 

Oscillatory divergence 

Fig. 5g 


response in Eq. (31), while discontinuous, is reasonably 
well approximated by an exponential decay. It is then 
convenient to describe the system’s response characteristic 
in tenns of a time constant of an equivalent RC circuit. A 
value of r near unity corresponds to having the system 
time constant larger than the conical-scan period When 
r is much smaller than unity, the response in Eq. (31) is 
very discontinuous but a proper interpretation still allows 
a time constant to be defined. To determine the lime 
constant, we shall first define Icq as the number of scans 
necessary to reduce <f>i in Eq. (31) by 1/e. Thus, noting 
that A: - 1 in Eq. (31) is the number of scans to obtain 
we have from Eq. (31) the equation to solve for /cq 


ing the low edges of the "‘staircase.” For small values of r 
the staircase steps are large, and what could be considered 
the average decay is »F/2 to the right. When r is small 
the steps are small, and the actual decay is very close to 
the exponential decay. 

B. Performance in Noise 

In Eq. (27), and are independent random 
variables since depends only upon which 

are independent of because n(() is white. Thus from 
Eq. (27), 

Var[<j.* + ,]= r^Var[<^*]+Var[nl'>] (37) 

We wish to find the steady state variance of We do 
this by assuming that 

lim Var[<>fe]=o^ (38) 

>ac 

By taking Ac^oo in Eq. (38) and using Eq. (30), 




(32) By solving for we have 


Thus 


ol = hVj{\-r^) 


(40) 


and 


- 1=^0 log,'- 


^0= - '/log,'' 


(33) Solving for h from Eq. (28) and substituting in Eq. (40) 
gives 


(34) 






(41) 


Note that only fortuitously will be an integer. We now 
define the system time constant t as the amount of time 
corresponding to conical-scan periods 

T^koP=-P/\og^r (35) 


but since 1 -r^ = (l -(-r)(l — r), we have 


('-'•)ov 

{\+r)A^ 


(42) 


For future reference we solve Eq. (35) for r 

r = e\p{-P/r) (36) 

In Fig. 5(c) we show how the discrete steps match up 
with the exponential shape to define r for the case r = 0.7. 
We plot the exponential decay starting at (0, 1) and follow- 


In Subsection A, we commented that cases c and d of 
Table 1 converge absolutely at the same rate if the abso- 
lute values of r are equal. We now see there is a distinct 
disadvantage to having -l<r<0. From Eq. (42), since 
Oyy and A do not depend on r, we see that (I -/•)/(! -Hr) is 
monotone decreasing over -l<r< + l. Therefore, 
although case c with and case e with r = - rj will 

converge absolutely at the same rate, the variance of the 
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error in Eq. (42) is much larger for case e when Tq is not 
near zero. 

Using Eqs. (19) and (22) in Eq. (42) gives 


T^'(R) Vl+rj \BP] 


(43) 


We thus have the expression for the mean-square tracking 
error in the elevation axis. In Eq. (43) the first factor is a 
function of the source and operating temperatures and the 
antenna gain pattern. The second factor is a function only 
of r and thus, from Eq. (36), depends upon the ratio of the 
time constant and scan period. The last term shows that 
radiop ncr bandwidth and scan period decrease We 
could use Eq. (43) in this form but we choose to express 
the second factor of Eq, (43) in terms of t by using Eq. 
(36) and have 



TIME CONSTANT/PERIOD r/P 




OF 


Ts 


+g(li) 


g'(R) 





The second factor in Eq. (44) is well approximated by 
{P/2ry^^ for large t/P. Define the rate factor F equal to 
the ratio of the second factor in Eq. (44) to its approxima- 
tion as 


Fig. 6. Rata factor 


^ OP 


-^+g(R) 


g'(R) ^ 


(47) 


If r/P is not >2, Eq. (46) should be used with F taken 
from Fig. 6, 




l-e-^/r y/^ 1 


(45) 


Substitution of F {P for the second factor in Eq. 
(44) gives 


f OP 


+«(^) 




¥ 


(46) 


A \Aoi of F is given in Fig. 6 and it is seen that, except for 
very small t/P, F is a second-order correction to account 
for the lobing rate vs lime constant. 

It is anticipated that t will be equal to at least a few 
times P. Thus we may take F= 1 for t/P >2 and have 
the useful result 


We now have a^, which is the standard deviation (rms) 
value of the elevation angle error of the conical-scan 
tracker. The dimension of is the same unit used in 
calculating g'(P), i.e., normalized power gain per degree or 
per radian. 

The results in Eqs. (46) and (47) are very general and 
apply for any value of R. We may now optimize the 
tracker by finding the value of R which minimizes 
(holding everything else fixed). Since R enters Eq. (46) 
only in the first factor, we need only to minimize this 
factor. A special case of very great interest occurs when 
T^^Tqp- Wc may then ignore g(P) in Eq. (46) since it is 
< 1. Thus, is minimized when R is chosen where g\R) 
is maximized, i.e., where the beam pattern has the greatest 
slope. Remember, of course, this must be calculated from 
a plot of g(P) with power gain presented linearly and not 
in dB. 

In general, we may find R to minimize the first factor in 
Eq. (46) by differentiating it with respect to R and setting 


t 


7 


} 

( 


8 
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the result equal to zero. The lesult is the following equa- 
tion which must be solved for R 


I 




* OF 


+«(«) 


g'm 


(48) 


A graphical argument can show that \1 g{R) is monotonic 
decreasing for R>0, the solution of Eq. (48) gives an R 
greater than or equal to the value which maximizes 
When optimization of pointing accuracy is paramount 
one would choose R as above. However, there are other 
considerations, prime among them, the crossover loss 
g(/?). When the system is operating primarily for 
radiometric purposes the system will be increased 
by \/g{R), We shall shortly consider this tradeoff in detail 
for a Gaussian beam pattern. 

A second point is the conical-scan rate. When R is 
increased, the scan rate must be reduced so that mechani- 
cal stresses remain acceptable. According to Eq. (46), an 
increase in P has no first-order effect on a^. However, one 
effect we have not yet considered is system gain fluctua- 
tion. The slower we lobe the beam the more we are 
susceptible to gain fluctuation errors. The analog to a 
Dicke radiometer is clear, whereby we wish to “switch” 
(rotate) faster than the gain fluctuations. 

C. Performance Tradeoff 

Solution of Eq. (48) is probably only of academic 
interest since, in practice, a tradeoff of performance vs 
g{R) will have to be made. Probably the most useful thing 
we can do here is to calculate as a function of R, Tqp, 
r^, B and t and give the tradeoff vs g{R). When this is 
done, the solution of Eq. (48) can be taken directly from a 
graph, if desired. Use of an exact g{R) for the antenna in 
question would be ideal. However, the main lobe of most 
antennas is approximated quite well by a Gaussian gain 
pattern. We shall thus assume the normalized beam power 
gain pattern to be 


g()8) = exp[-^()8/»0^] (49) 


where IV is the antenna beamwidth between half-power 
points (with the same unit as R) and 

/i = 4 log, 2 = 2.773 (50) 

We now assume that rfP>2, so we may use Eq. (47). 
Thus, keeping only the magnitude of g'(/?). 
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NORMALIZED SCAN RADIUS R/W 

Fig. 7. Angular tracking error for radio source, gain fluctuations 
ignored 
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T 

' 5 


luR^Br 


( 51 ) 


Normalization of and use of Eq. (50) allow us to plot 
this relationship as 


iV 


= (0.180) 


T 
‘ s 

R/W 


(52) 


This is the major result for the standard deviation of the 
elevation error of the conical-scan tracker. The units of 
o^, R and IV must be the same. By symmetry, this result 
must also apply for the cross-elevation axis. In Fig. 7 we 
portray the rms angular error of Eq. (52) in terms of R/\V 
and Tqp/T^. Note carefully that i^/lT=0.5 corresponds 
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to the half-power point of the antenna beam. The system 
crossover less is the system gain on boresight vs the 
antenna beam peak gain. From Eqs. (49) and (50) it is 
easy to show that the crossover loss in dB for a Gaussian 
beam is 

g(R) in dB= \2M{R/WY (53) 

We show this degradation in Fig. 8 so that the tradeoff of 
tracking error in Fig. 7 and crossover loss can be seen. 

In Fig. 7 i t is easy to see the values of R/W which 
minimize o^^Br /W for fixed values of Tq^/T^. When 
simple calculus shows that g\R) is maximized 
and hence is minimized for 

l/V^ =0.425 (54) 


The derivation above of has not included system 
gain fluctuations. Just as in Dicke-type radiometers, if the 
switching rate (here we have rotation rate) is not faster 
than the gain fluctuations, the performance will be limited 
by gain fluctuations. It is anticipated that scan periods of 
10-100 seconds arc too Mow to overcome the gain fluctua- 
tions. Thus, the result of Eq. (52) is optimistic, especially 
for large values of the s;an period P. In Section V of this 
analysis we shall treat gain fluctuations in the systems of 
interest. 

It should be noted that as we have derived it applies 
to each antenna angle axis. The errors in these axes can be 
assumed to be indepimdent so the total instantaneous 
angle tracking error (radially) is 


c =- ^ 9 ^ 4- 


From Fig. 7 this approximation for the minimum point is 
seen to hold very well down to at least Tqp/T^^ 10. 



0 0.2 0.4 0.6 0.8 ].0 


NORMALIZED SCAN RADIUS R/\V 

Fig. 8. CroMover Iom 


Because /i(r) is a Gaussian process, 0 and 6 are Gaussi:.n 
since the system is linear. Thus the radial tracking error c 
is a Rayleigh random variabK; with probability density 
function 





lo 


€>0 

€<0 


(56) 


The average radial error between boresight and target 
thus is 


€ = £[«]= f“y(e)d'€ = (^/2)'/\=1.25a^ (57) 

•'0 

D. Operational Considerations 

To operate the system considered here, various choices 
of T (or equivalently r) must be selectable. From Eqs. (35) 
and (28) we have 


r-^-P/\og^(\-hA) (58) 

We also have A in Eq. (19), 

A^CBT^-{R)P/2 (59) 

Therefore, we see that r is determined by P, h, C, B, Tg 
and R (via g\R)Y Even if P, h, C B and R were held 
fixed, we see that the system time constant is a strong 
function of In order to maintain constant t for 


;>. PH0DUCIBILITV op the 
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different radio sources, the gain parameter h must vary 
inversely with Ty. In addition, if the scan period F, the 
radiometer gain C or bandwidth B, or the scan radius R 
change, then k must vary so that t will stay fixed. In most 
tracking sy steins there exists an AGC loop to maintain 
constant track ng bandwidth over a wide range of signal 
levels, kl ! H ent, there is no AGC loop for ladiometric 
work, >0 h vill have to be set according to a calibration 
schema wh h takes each of the above parameters into 
account. Solving Eqs. (58) and (59) to give h for fixed t 
we have 

A *2(1 - e-^^^)/CBTsg'{R)P (60) 


pattern to be well approximated by a Gaussian shape as 
in Eq. (49). The normalized antenna voltage gain is thus, 
by square-rooting (Eq. 49), 

= (62) 

where 

|u' = fi/2= 1.387 

The normalized voltage gain vs time is, from Eqs. (4) and 
(62), 

e(0 = e[/?(0] 


Within the r mputer program, calculation of h is readily 
made when \ t, C, B, and R (assuming the program 
has the g'{R; function) are fed in. The value of C can be 
meas\ircd (or adjusted to a prescribed value) by observing 
the radiometer output when observing the system ambient 
load or cold ;.ky. The value of B is usually close to some 
nominal value for a particular radiometer. The value of P 
may be set up as an adjustable parameter or it may be 
fixed in the program. Thus, t and are left as variable 
parame^ ;rs for system operation. A direct calibration to 
measure A of Eq. (19) for a particular source may be a 
convenient technique and has the advantage of checking 
all syste ,i parameters ir: one measurement. This could be 
done by observing the integrator outputs, Eqs. (15) and 
(16), wht 1 the system pointing is running on predicts plus 
lobing alone, with the offs^*' loop open and deliberate 
Known errors in 9 and ((> _„iSting. One problem with this 
method, howtwer, is in knowing 9 and <f> exactly if the 
exact source position is not well known. 


+ 2y?(tfcosw^/ + <f>sinw^/)] I (63) 

We then see that the signal voltage received when the 
antenna is conically scanned is 

r(/) = i;(/)5(0 (64) 

At the point where we define signal power we also 

define the system operating temperature ^OP which gives 

N^^kTop watts/Hz (65) 

on Zq ohms, where k is Boltzmann’s constant. 


IV. Spacecraft Tracking 

In this section we will derive the performance of the 
conical-scan trad er when i* is used on spacecraft. We 
shall ignore gam fluctuation 'Jects and defer their con- 
sideration to Section V. "'^'hen the antenna’s peak gain is 
pointed at the spacet kt, we receive a carrier power of P, 
watts in the antenna. This occurs as a signal voltage 

cos (wq^ -♦■»//) (61) 


on an impedance of Zq ohms, at some reference point, 
'fhe quantity \p is the ohase of the received signal. The 
lobing procedu’^e in Section II produces the instantaneous 
angle in Eq ;t) by which the antenna’s peak gain is 
away from he spacecraft. We assume the antenna gain 


Without loss of generality we may take Zq — I ohm and 
then the received signal plus noise is 

r{t) = v{t)yj2P^ cos ( wo / + i |') + m (/) ( 66 ) 

where n{t) has the two-sided power spectral density 

S„U)^No/2 (67) 


A. Manual Gain Control 

We shall first consider conical-scan tracking when the 
receiver AGC is not used and the receiver gain is fixed. 
We treat this case for two reasons: (1) it can be solved 
exactly and provides a, jumping-off point for the 
approximate analysis of the AGC case and (2) when the 
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AGC is faster than the scan rate the manual I’ain control 
(MGC) result will apply. 

We must detect the amplitude of the signal in Eq. (66) 
so that we may correlate with sinw^/ and coso^/ to give 
pointing information. Coherent detection is accomplished 
in a phase detector used for AGC detection by means of a 
reference signal generat^^d by the phase-locked loop 
(PLL), The PLL tracks the phase of the sgnal and 
gt.ierates a reference signal 

0 = ^ cos(w„^ + ) (68) 


where 

/•/* 

I (75) 

•'0 

= a: f d/ (76) 

•'0 

For simplicity, these integrations are tacitly assuming a 
time scale starting at r « 0 at the beginning of each scan. 
From Eqs. (63) and (75), 


where \p is the PLL’s estimate of tp. We ignore he erodyn- 
ing and equivalently assume the PLL to be generating its 
reference in the RF range. The phase detector forms 

>(0-e„K0K0 (69) 

Ignoring double-frequency terms we have 


4»exp 





fji 

IT' 


cos «„/ + <).* sin V) 




(77) 


/(0 = u(/)V/’,cos(^-«^ ) + «'(') (70) 

where n'(0 is again a white process with power spectral 
density (Ref. 5) 

5..(/) = ATo/2 (71) 

If the PLL loop SNR is above about 6 dB. ^ is close 
enough to V' that cos(V'-’V')^=^L We assume this to be the 
case here. If the PLL is operating nearer than 6 dB to 
threshold this approximation will be optimistic. We thus 
have 

y(i)^v(t)\/^ +n'(t) (72) 


The exponential under the integral can be re';»'ritten as 

exp[-X*cos(«„/-A*)] (78) 

where 


+ /IV^ (79) 

and 

•'<-“"''(1) w 


Because of receiver gains the actual signal presented at 
the signal strength point in Fig. 2 wil! be different by 
some gain K over Eq, (72). Including this in Eq. (72) we 
have the signal strength jc(/) 


The Jacobi-^Anger formula is (Ref. 6) 

exp(Xcos ft) = /o(X) + 2 S /n(4)cos nft (81) 


x(t)^K\v(t)\f^ +n'(0] 


(73) 


where J„(i) is the modified Bessel function of the first 
kind. We may expand Eq. (78) using Eq. (81) and obtain 


We now use Eqs. (15) and (16) to develop the offsets for 
the elevation and cross-elevation axes. We start with ele- 
vation and have 


Ar»=-/,[Arv'/’.4< + ^l'’] (74) 


exp[-X*cos («„r-\*)J 


“7o(-X*) + 2 2 7,(-4*)cos [«(a'„/-A*)] 


n- 1 


(82) 


I 

i 

1 
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Use of this in Eq. (77) gives on’y the following te’tr. since, to W. Therefor the argument of /| in Eq. (87) will be 

for # 1, the sinusoids in F.q. (82) are orthogonal to much smaller than unity and we may use the first term of 


sinw^f: 






rP 

X I sin u>„( cos(w^/ - \)dt 
•'0 


(83) 


the power series expansion of /i(x), 
/,(x)asx/2 

Using this fiq. (8"0 becomes 


( 88 ) 


PiiR 






(89) 


4"^’2,( -$(i)si'iX*exp[ - <#>;[) 

L ^ 

From Eq. (80) we see 

cl 


Also 


Thus 


/,(-x)=-/,(x) 


(84) 


(85) 


( 86 ) 


;=--=xp| 




+ <t>l 


(87) 


The result in Eq. (87) is exact, but is more precise than 
needed.* We may simplify this result with essentially no 
loss in accuracy. Note that 






is the radial error. This error will generally be much 
smaller than W. Also, R will be less than or comparable 


Further, we assume and \ <(>^| to be «:/?. We then may 
neglect and in the exponential and we have 

4=-<>*~exp[-M'(^/H'r] (90) 


This last approximation will be valid unless R is chosen 
very small and and are comparable to R, 

From Eqs. (74) and (90) we liave 

(91) 

where 


KPy.Ry[p, 

A ^ (92) 


We dropped the minus sign in Eq. (90), so the correction 
is in the right direction. In Eq. (91) we have the same 
equation as Eq. (17). We need only find the variance of 
and we shall be able to directly use Subsections A and 
B of Section III to find performance. 

We know has zero mean because n\t) is white. 
Also in the Appendix we substitute Kn\t) for «(/) and 
obtain 


•'0 •'o 


' In Section III, linearization was assumed first and then a Gaussian 
beam. Here a Gaussian beam was assumed first and now linearization 
is to be done. Either sequence could be used in either place— the result 
is the same. 


Since n'(() has the power spectral density in Eq. (71) we 
know 




(94) 
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Use of Eq. (94) in Eq, (9 ^) gives 


ol^PK^NJA 


As in Section III we argue that F«=:l for practical applica- 
tions. Also by use of Eq. (65) we obtain 


(95) 


We now can use Subsection A of Section 111 in its 
entirety for the spacecraft case. We may also use Subsec- 
tion B through Eq. (42) and the paragraph immediately 
following it. Upon substitution of Eqs. (92) and (95) in 
Eq. (42) we have 


exp[ V^'{R/Wf] 


(98) 


y/PnR\fK 


(96) 


Normalization and putting in the values of /yi and ft' give 

(99) 



Expressing r as in Eq. (36) and defining the rate factor F 
as in Eq. (45) we have 




(97) 


By symmetry this result also applies to the cross-elevation 
axis as well. This result is presented in Fig. 9. Boltzmann’s 
constant k has been included in the normalization be- 
cause it is felt that ^/FJkT^ is most easily computed 
using dB, whereby A is -198.6 dBmW- K- Hz. Note 
that Pjr/kTQp is dimensionless bec uise 1/r has the di- 



Flg. 9. Angular tracking error for apacacraft, gain fluctuations 
Ignored 


Fig. 10. Tracking error vs crossover losi, gain fluctuations 
Ignored 
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KPy.R 


( 101 ) 


in Fig. 8. F^r ease in comparing vs ^/?) we show their 
tradeoff in Fig. 10, which is a cross-plot of Figs. 8 and 9. 

If all parameters are held constant except R, simple 
calculus shows that the minimum in Eq. (83) and Fig. 9 
occurs when /?-(0.601)W^. If pointing accuracy were 
paramount then this choice of R would be optimum 
unless the crossover loss of 4.3 dB given by Fig. 8 causes 
the phase-locked loop to approach threshold. From Fig. 
10, it can be seen that permitting to be only twice its 
minimum attainable value allows ^/?) to be reduced from 
4.3 dB to 0.45 dB, a sizable reduction. 

The reader must be cautioned that gain fluctuation 
effects have not been included in this calculation but will 
be considered in Section V. The present result is sufficient 
to characterize performance as long as the conical-scan 
rate is faster than the gain fluctuations. 

B. Automatic Gain Control 

To maintain signal levels in the PLL receiver near 
nominal values over a wide range of input signal levels, 
AGC is used. A secondary benefit is that of maintaining 
approximately constant closed loop bandwidth when per- 
forming a tracking function such as monopulse pointing 
or linear polarization tracking. This advantage also occurs 
in the conical-scan tracker. However, it is a mixed bless- 
ing since the choice of AGC speed will often not be made 
to optimize pointing but to optimize PLL locking per- 
formance. If the AGC is very slow the error signal we 
develop by conical-scan tracking will be attenuated and 
phase-shifted by the AGC circuitry. Let us begin by 
considering the AGC characteristic with respect to the 
steady component of signal level (ignoring temporarily the 
small sinusoidal error signal). It is known that the 
dynamic AGC voltage of the PLL receivers used is very 
close to - 1.0 V over a signal range of better than 60 dB. 
Thus it is clear that the receiver gain must vary inversely 
as the signal voltage. Thus, when using AGC the K in Eq. 
(73) will actually be 


( 100 ) 


where K' is a fixed quantity and is, of course, 
proportional to signal voltage. If the AGC response is 
much faster than the conical-scan rate so that the 
AGC gain at frequency is the same as for the DC 
steady stale signal component, we argue that K^K'/yfp^ 
also for the error signal and Eq. (92) becomes 


2W^ 

Also, from Eqs. (86) and (93), we now see that the system 
closed loop time constant theoretically does not depend 
upon the signal level P,. This represents a distinct advan- 
tage over the conical-scan tracking of radio sources in 
Section III, where t does depend upon the source tem- 
perature T,. 

The performance of the comcal-scan tracker under the 
assumption that the AGC is much faster than is just 
that given in Eq. (99), since using Eq. (100) in Eq. (95) 
and by use of Eqs. (101) and (42) again gives Eq. (96), and 
everything down to Eq. (99) follows. 

In practice, however, the AGC may not satisfy the 
criterion of being much faster than <o^. We now consider 
this case. For small signal aii.plitude variations (appli- 
cable to the error signal) the AGC circuit acts as a single 
pole filter with transfer function (Ref. 7): 

where b is dc gain and is the filter response time. Note 
that is not the so-called “AGC filter time constant” 
which is «380, ?»34, or j5s:4s in DSN Block III receivers. 
These numbers represent the time constant of an RC 
operational amplifier circuit approximating an integrator 
in the AGC loop. Since these lime constants are in the 
feedback loop, the overall response time is determined 
by the open loop gain, which is a slowly varying function 
of signal level. To find t^, it is a simple matter to put a 
small step change in a test transmitter signal level and 
observe the AGC response time to 1/e. For example, a 
filter time constant of 380 s generally gives a value for 
which is on the order of 5 s. 

If ^ 1 /w^, the effect of the AGC is nil since the error 
signal is affected the same as the steady signal component. 
This is the case we treated above. However, if increases 
toward I the AGC causes attenuation and phase shift 
of the error signal. The apparent error signal amplitude 
will be diminished and the tracking response will thus be 
slower. This arises due directly to the AGC attenuation 
and indirectly to the AGC phase shift because the error 
signal will be out of phase with the reference signals in the 
correlation procedure. One way to resolve the matter is to 
apply analog compensation to the AGC voltage before 
performing the correlations. A suitable filter would be the 
following: 
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Wcomp(w)” 1 


(103) 

which cancels the pole in Eq. (102). Practical implementa- 
tion dictates that £q. (103) has, in addition, a second- 
order high-frequency roll-off determined by the analog-to- 
digital sampling rate. With the compensation of Eq. (103), 
performance of the conical-scan tracker would be identi- 
cal to that of the fast AGC case above, i.e., the MGC 
result in Eq. (99). The compensation in Eq. (103) can also 
be accomplished digitally in the computer after sampling. 

A second way to remove the smoothing effect of the 
AGC is to apply a gain to the error signal equal to that 
loss caused by the AGC and apply a phase shift to the 
correlating reference signals equal to that caused by the 
AGC. If T is several times P so the error signal changes 
slowly (thus looking approximately like a sine wave) the 
system response will be equivalent to that obtained by 
using the analog or digital compensation discussed above. 
The reason is that if the error signal looks like a sine wave 
with slowly varying envelope the system will be looking at 
the noise spectrum only for a small bandwidth fttl/r 
around Over this small bandwidth the AGC transfer 
function is approximately constant. Hence, the error 
signal looks like a narrowband signal plus white noise, 
and the effect of the AGC is equivalently a gain and 
phase shift at When t and P are of comparable 

magnitude this argument does not hold because the error 
signal no longer is narrowband, and the analysis is very 
difficult and will not be carried out. Logic indicates, 
however, that the system will still work and it is felt that 
system performance will not significantly change from the 
T<P case even up to t“P. Since it is anticipated that 
T<P, it is felt that an exact analysis for arbitrary t and P 
is not really needed. 



NORMALIZED TIME C NSTAr."^ P * w r. ' 2tt 
A ff) A 


Ftg. 11. AQC gain magnNuda 


For ease of 'calculation, the attenuation and phase shift 
of the error signal due to the AGC arc presented in Figs. 
11 and 12. Gain magnitude is 

-^^-1/Vi+(2wt^/P)^ (104) 

and phase shift is 

Angle[W(«„)]- -tan-'(2^r,//’) (105) 

From Figs. 1 1 and 12 it is seen that t^/F*0.05 reduces 
the gain by only 0.95 and gives only - 18 deg of phase 
shift. Remembering that phase error causes the output 
amplitude of the correlation procedure to drop as the 
cosine of the . . ror, the total drop in loop gain (vs the fast 
AGC case) is 0,95 cos (-18 deg) *0.90, a small change. 
Thus, as long as t^/F<0.05 we probably do not need to 
worry about compensation. Note that Figs. 1 1 and 12 and 
the argument above apply equally well to the specification 
of the “fast” RC filter when tracking radio sources (see 
Section III) and was the argument used there to require 
that T/jc/F< 0.05. It should be observed that a phase 
error of 30 deg or more gives substantial “cross-talk” 
between elevation and cross-elevation axes, which results 
in a deterioration of performance and can produce an 
unstable system as the phasing error approaches 90 deg. 

In conclusion, we find that the performance under 
AGC is the same as for MGC, namely Eq. (99), as long as 
we compensate for the AGC filtering action when nec- 
essary. 



NORMALIZED TIME CONSTANT P w 

A m A 


Fig. 12. AQC ph«M shift 
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V. Gain Fluctuation Effects 

In Sections 111 and IV, the analyses of tracking per- 
formance were made under the assumption that fluctua- 
tion of the system gain was of no consequence. In a 
practical case, this would be justified only if the scan rate 
were much faster than any significant spectral component 
of the gain fluctuation. In this section we shall consider 
the gain fluctuation problem in detail. 

The gain of all receiving systems varies with time. 
Mechanical flexure due to vibration, varying stresses, etc., 
and thermal variations are the well-known primary causes. 
Good receiver design usually can reduce the problem to 
an insignificant level everywhere except in active devices 
at RF frequencies. For DSN receiving systems, the masers 
are thought to be the principal contributors to gain (and 
also phase) instabiliues. The gain stability is excellent by 
most standards and is satisfactory for virtually all com- 
munication purposes. However, for radiometry and the 
present tracking applications, the gain fluctuation of the 
system will be seen to give the fundamental limit for the 
tracking accuracy available. 

We now define the «eceiver gain as the gain from the 
feed horn down to the point where we convert to digital 
format as 

e(0-(2o[i+*'='(0] («06) 


where f{t) is the gain fluctuation. We really do not need 
to specify Qq very carefully, as it will drop out of the 
analysis. All we need is to realize that the system gain 
varies as 1 + F{t), For good receivers, we have 

\F{t)\<^\ (107) 


Thus the voltage gain clearly varies as 
I ^(0 

Vl+F(/) ^1 + ^ (i08) 

We shall model F(t) as a stationary random process and 
we assume that we have knowledge of the statistics (i.e., 
spectrum) of F(/). 

A. Tracking of Radio Sources 

We start with Eq. (14) and multiply it by our time- 
varying gain [1 +F(01 to give x'(/): 


x'(0-[l + ^(0]^(0 (109) 

-CB[l + F(l)][Tof + T,g{R)] 

+ CB[ l + F(l)]TsgXR)[0 cos«„/ + <»sin«^] 

+ [l+F(f)]n(0 (110) 

Since |F(/)1<1, we may neglect F{f) with respect to unity 
in the second and third terms of Eq. (110). We may do 
this because in the second term F(i) causes a very small 
random variation in closed loop bi.nd width and in the 
third term F{t) causes a very small random variation in 
the input noise level. Neither of these are first-order 
effects. In the first term of Eq. (110), /*X0 <^oes cause a 
first-order effect, as we shall see. Upon neglecting F)(/) as 
above we find from Eq. (110) that 

x'(0-af(0+x(0 (111) 

where 

A-CB[r«, + 7-sg(F)] (112) 

We now form and Ajf^ as in Eqs. (15) and (16) by 
substituting x'(0 for jc(/) therein 

X'{t)smw„tdt (113) 

•'(A - \)P 

a / x'{t)cosu)^t di 

^(k-\)F 

We carefully specified the integration limits in Eqs. (113) 
and (114) instead of starting /*0 at the beginning of each 
scan as in Section HI because we will need to consider 
correlation between the values of F{t) in different scans. 
We now will consider just the elevation axis and appeal to 
symmetry later. By putting Eq. (Ill) in Eq. (113) and 
integrating as in Section 111, we get Eq. (17) plus an extra 
term which is due to the first term in Eq. (Ill) 

("5) 

where 

F(/)sino)jdf (116) 

•'(* i)/» 
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( 122 ) 


A is as in Eq. (19) and is as in Eq. (20). Starting with 
Eq. (17), wc follow down to Eq. (30), keeping the extra 
term in Eq. (115) added on. We then have in analog 
to Eq. (27) 




■£ 2 2 

pmO 


pmO 1-0 

(12J) 


where r and arc specified in Eqs. (28-30). 

Because and are not independent when j ^ k, 
wc shall have to consider the actual solution of Eq. (117) 
rather than just evaluate variances as independence al- 
lowed us in Section 111. The exact solution of Eq. (117), 
assunrdng wc start at /»0 with is 

A- 1 

2 + (118) 
p^\ 

Note that the summation term is random while the first 
term is just Eq. (31). Thus wc sec the response time 
characteristic is the same as in Subsection A of Section 
III. Wc arc primarily concerned with the steady state 
variance of so wc take k very large. Since |r|< 1 ihc 
first term of Eq. (118) drops out. Wc make the substi .u- 
tion i^k-p- \ in Eq. (118) and have 

1-0 

Since and are stationary, we may drop k- 1 from 
the two subscripts in Eq. (119) and reverse the sign of / 
without any change in the statistics of *hcn have 

* 2 (> 20 ) 

1-0 


Now from Eqs. (28-30) and the first sentence of Subsec- 
tion B of Section 111, 

{*‘”j '-f (124) 

and 

(125) 

Thus 

<t^”/i^o^2 '■^'+ 2 r‘^PE[m<:‘W;^] (126) 

1 — 0 j — 0 ^ — 0 

The first summation is 1/(1 -r^) so the first term of Eq. 
(126) is h^al/(\ - which is just in Eq. (40). Thus 
the first term of Eq. (126) is just the steady stale variance 
of when F(/)*0. We then write 

+ (127) 

where is the second term of Eq (126) and is the 
variance due to gain fluctuations, and is as calculated 
in Section III (or Sectu>n IV when discussing spacecraft 
tracking later in this section). 

We now need to find a^. From Eqs. (126) and (116), 


Wc desire the variance of as Ac— >oo so if we take k-^oo 

in Eq. (120) wc have 


2 2 £[£(/, )£(/2)]sinw„r 


;-0 />-0 


1-0 

The steady state variance of is then 


[ 121 ) 


X sin u„t 2 dr, dr 2 


(128) 


The expectation in Eq. (128) is the autocorrelation of 
F{t) 


o^-Var[<^.^]-£[<})^^J 


E\F(t,)^\t2)]’^R^(t,-t2) (129) 
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which may be expressed in terms of its power spectral 
density Sp{f) as 

/“ (130) 


The first summation is 


2 ) “ ( 1 “ ) * 


(138) 


1-0 


where 


and the second one is 


y-V-1 


(131) 


2 

p^Q 


-1 


(139) 


Substitution of Eqs. (129) and ^130) in Eq. (128) ami 
interchanging integrals gives 


The product of the summations can be maniniiated to 
give (l+r^-2r cos 2^fPy \ We can thus rew*i*e Eq. 
(137) as 


1 — 0 />-0 


(132) 


where 


[«i,-(2^/)'] (I + 


r M/)yp(/)yr(/W (133) 

- 00 

where * denotes complex conjugation, and 


— d/ 
~.2n/P) 

(140) 


•'(/>- 1)/* 


Evaluation of y.(/) via much algebra gives 


(134) 


Upon substituting h~(l - r\J a from Eq. (28), A from Eq. 
(1 12) and A from Eq. (19) and using 

f„~uJ2^~\/P (141) 

where f„ is the conical-scan rate in Hz, we have 


Y,(/)' 


- 1 ) 

«^-(2^/)' 


. T, 


(135) 


OP 


+K(«) 


K(fi) 


(142) 


Substitution of Eq. (135) into Eq. (133) gives 


where 




S;: (Z)sin^ f^df 


(136) 


r~ f" Sf (fm/)df (143) 


Substitution of this into Eq. (132) and rearrangement and the weighting function H^(/) is 
gives 




[wi.- (2^0^] 

X f ) 2 )^df (137) 


W(fy 




^ 1 




(144) 


1-0 


/» — 0 


[l-t-r^-2rcos(2w///„)] 
By use of Eqs. (127) and (47) we obtain 
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If 



a« 



+«(«) 


g’iH) 



We note for r near unity that the cosine function in the 
denominator of Eq. (144) is the principal factor determin- 
(145) ing the area of W(f), We thus expand the numerator of 
Eq. (144) around /-/„ 


'Dus is the principal general result for tracking of radio 
sources. 



It must be realized that T from Eq. (143) depends on r 
via W{f) in Eq. (144) and hence upon t, since from Eq. 
(36) 

(146) 


(148) 

Also, the first term in the denominator of Eq. (144) can he 
similarly expanded around /»/„ as 


We now wish to examine in detail how F depends upon r 
(via r), and S^(/). In Fig. 13 we plot from Eq. 
(144) as a function of / for v'*-ious values of r. It can be 
shown that the W{f) for any r-|rol is smaller than that 
for /■•“jrQj. Thus, we are further justified in taking 
r6(0,l) as wc argued initially in Section III-B. In general, 
we must calculate F from Eq. (143). However, in the next 
section we shall show that in most cases of interest, we 
can approximate F quite well by a simple expression. 



(149) 


Using Eqs. (148) and (14';, in Eq. (144) wc find W( f) near 
/„ to be given approximately 


W(f)^ 


(1-^f 

l+r^-2rcos(2ir///„) 


(150) 


B. Narrowband Approximation 

Wc sec in Fig. 13 that for values of r near unity that the 
area of W'(/) becomes concentrated near /*/„. We know 
that in almost all cases of interest wc will have F/t<1 so 
that r is near unity. Wc thus realize that in these cases wc 
can approximate W(f) with delta functions at /“/„ and 
/• if Sf{f) is reasonably smooth near This 

smoothness has been observed, so we are justified in 
taking >F(/) as a pair of delta functions, Wc now need to 
find the area near /“/^ in W{f) so wc can approximate 
Wif) as 

nf) - -U (147) 



0 0.5 KO ) 5 2.0 2.5 3.0 


NORMALIZED FREQUENCY 1/1^ 

Fig. 13. Walghtlng function 

20 


Now since Wq is to be the area of W(f) near wc 

can find Wq by integrating over one period of Eq. (150) as 

W'o- (151) 

This integration can be carried out and after much manip 
ulation gives 

(152) 

From Section III-B, wc know (I -r)/(l +r) is well 
approximated by P/2r. Using this and the obvious rela- 
tion F- l//^,,wc have 

W^o-l/2r (153) 

Thus, from Eq. (147) 

[«(/+/..) + «(/-/.)] (154) 

and using this in Eq. (143) wc have 
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r-^[s^(-/J+M+/«)] (»55) 

Sfif) is an even function, however, so Sgr{-f^^Sf 
(+/«) we obtain the simple result 

r = 5f(/„)A (156) 

We refer to the approximations made in this section as a 
narrowband approximation since by restricting t to be 
large the system is responding only to the gain fluctuation 
noise in a narrow band around Substitution of Eq. 
(155) in Eq» (145) gives 


to evaluate o. It is interesting to note that a varies 
inversely with t, for arbitrarily large values of t. In 
addition, we see that gain fluctuations p>revent one from 
achieving arbitrarily small a by making B very large. We 
see, in fact, that if 1/B is comparable to 5^(/,J we can 

decrease a only by an additional factor of even if 
Since S^if) is almost always monotonically de- 
creasing, we can also see mathematically the decrease in o 
which can be obtained by increasmg/^. It is necessary, of 
course, to know Sjr(/) before c can be calculated. 

We have modeled the beam snape as Gaussian so the 
crossover loss was found in Eq. ( 53) to be 


T 

r * / 


OP 


5 




\/2 


r »/2 


(i^7) 


g{R) in dB- \2,04{R/Wf (161) 


and was plotted in Fig. 8. In Fig. 15. we present o in Eq. 
(160) cross-plotted against the crossover loss g(R). In this 
plot it is easy tc see the tradeoff of a vs g(/?). 


We now take g( p) to be Gaussian as in Eq. (49) 


^(^) = exp[-M(iS/»n'] (158) 


Substitution in Eq. (157) gives 


1 + 

' S ' 


»/2 





(159) 


Use of Eq. (50) for and writing this in normalized 
fashion gives 


or 


1/2 


m 




T/2 


= (0.180)- 


rJw 


(160) 


This noimalized function is the same as that in Eq. (52) 
except for the presence of Sp(f^). so Fig. 7 can be used 
for evaluation. We present this result in Fig. 14, We thus 
have the rms error in the elevation axis as a. By symmetry 
the error in the cross-elevation axis is the same. The 
simplicity of this result is impressive. We see that when 
T»P and S^{f) IS smooth, only is needed in order 
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C. Tracking of Spacecraft 

We assume as in Section IV that the receiver is either in 
MGC or compensation has been made for the effects of 
AGC upon the error signal. We saw in Eq. (108) that the 
system voltage gain varies as [l+F(;)/2] so we take the 
signal strength jc(/) in Eq. (73) and multiply by this factor 
to obtain 

x\t)=[\+F(t)/l]x{t) (162) 

= a:[h-f(0/2]V^i^(0 

+ [l+f'(0/2] Kn(t) (163) 

We can ignore F{t) in the second term of Eq. (163) since it 
gives only a ver>' small variation to the effect of n\t). In 
the first term of Eq. (163) we can determine an average 
value of o(0 as far as F{() is concerned. Looking at Eq. 
(63) and realizing that usually then 

average[u(/)]a£exp[ -ju'(/?/H^)^ j (lb4) 

(If B and are not <r/?, we will usually be able to assume 


the average is ;^1.) Using the above approximations, we 
have 



x'(0=af(0+^(0 

(165) 

where 

A = ^^^exp[-M'( /?/»')'] 

(166) 

We now form and A^"^ as in Eqs. (113) and (114), 

using the argument after those equations to write the 
integral limits as [(A - \ )P,kP\. We obtain for the eleva- 
tion case 



(167) 

where 

fkp 

•'(A-1)/* 

(168) 


n\t)s\Tu^jdt 

7(a-1)/> 

(169) 


—h\ f F(/)sinoj^f df 

7(a- i)p 

(170) 

In Section IV, using Eqs. (77)-(90), ia evaluated so we 
can rewrite Eq. (167) as 



(171) 

where, as 

in Eq. (92), 



KP^R\Ps r ,1 

2IV= “'’t '■'<''''"'>1 

(172) 


The are independent random variables with variance 
aj given in Eq. (95). Equation (171) is now in exactly the 
same form as Eq. (1 15), with in Eq. (170) identical to 
Eq. (116). The only difference is that A is given by Eq. 
(172) rather than Eq. (19) a. id o\ is given by Eq. (95) 
rather than Eq. (22). We may thus follow our previous 
work from Eq. (1 15) on down to Eq. (140). 
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1/2 


We now substitute h^(\-r)/A from Eq. (28), A from 
Eq. (166), and A from Eq. (172) and using Eq. (141) wc 
have 


Use of this in Eq. (174) gives 






4j[i" R 


(173) 


fiRy/l'T 


S,{fJ 


( 176 ) 


where F is as in Eq. (143). From Eq. (127), where is 
taken from Eq. (98), we have the rms tracking error given 
by 


We now put o in normalized form and use Eq. (50) for jit 
to obtain 


a — 




1/2 

(174) 

a 

T 

liR 

P^lr 4 


IF 

SfLL) j 



= (0.180) 


As in the radio source case, the general case involves the 
evaluation of Eq. (143) for F. However, since r will be 
near 1, we can invoke the narrowband approximation of 
Section V-B and have, from Eq. (156), 


(177) 


This result is presented in Fig. 16. We also present this 
result in Fig. 17 cross-plotted against the crossover loss 
g(R) in Eq.(161). 


r = S,iU/r (175) 



0 0.2 0.4 0.6 0.8 1.0 


NORMA'.IZeO SCAN RADIUS RAV 

Fig. 16. $pacacraft angular tracking arror 


It is of interest to compare the rms angle tracking errors 
for the radio source and spacecraft cases when the system 
signal-to-noise ratio becomes very large and only the gain 
fluctuations limit performance. In Eq. (160), if we take B 
and to be infinite, we get the same resu as when we 
take in Eq. (177). This is not a coincidence since 
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«ve should get the same result when the only limitation is 
the system gain fluctuation. 

VI. Experimental 'mplementatlon at the 
Goldstone 64-m>diameter Antenna 

In the period 1972-1974, an experimental implemcnta' 
tion of he conical-scan scheme was made on the 
NASA/JPL 64-m-diameter antenna Mars Station (DSS- 
14). The system worked well for tracking of both 
spacecraft and radio sources. The system has been used 
extensively in an operational quasi-manual mode for 
tracking since the initial installation. A fully automatic 
operational conical-scan program is now being developed 
for regular use in the Deep Space Network. 

A. Conical-Scan Hardware 

The conical-scan system is easily broken down into 
hardware and software blocks. The hardware diagram is 
shown in Fig. 18. The system permits any of five signals to 
be used as conical-scan ‘‘signal strengih’" data. The 
broadband square law detector (Ref. 8) can be used with 
any of the RF front-ends on the antenna. The four 
receivers at DSS-14 each have their AGC voltage brought 
to the Multi-Channel- Analog-Signal-Conditioner 
(MCASC). The MCASC consists of a separate DC ampli- 
fier, lowpass filtering and level shifting for each channel. 
The purpose is to bring each of the input lines into a 
range compatible with the analog-to-digilal converter 
(ADC). The ADC has 12 conversion bits including sign, 
and it samples 10 times per second. So that there is less 
than 4% decrease in SNR due to the finite sampling rale 
of the ADC, the time constant of the MCASC must then 
be O.l s and it is set accordingly (Ref. 9). It is also 
desirable that quantization noise not be a limiting factor. 


Thus there must be enough gain prior to the ADC, and 
the ADC must have enough levels that the voltage into 
the ADC randomly crosses several levels due to noise 
alone. A 12-bit ADC was found to be adequate for a wide 
range of signal types. 

B. Conical-Scan Software 

The computer used for the Antenna Pointing System 
(APS) at DSS-14 is a Scientific Data Systems SDS 910. 
The conical-scan program was packed into approximately 
2000 24-bit locations which were available after the 
primary APS program was placed into core. The conical- 
scan program is an option which is selectable by an 
operator by using a breakpoint switch. 

A block diagram of the computer calculation is shown 
in Fig. 19. The program was configured so it could do a 
conical scan in either hour angle /declination (HA/DEC) 
or azimuth /elevation (AZ/EL) coordinates. The 
appropriate secant correction is required in either case so 
as to get a circular scan as seen as a projection on the 
celestial sphere. The operation is easy to follow. Consider 
an AZ/EL scan. The scan frequency in radians/ second is 
w^ = 27t/P. The AZ scan /^sinw^t is corrected by the 
secant of EL and is fed to the antenna pointing function 
along with the EL scan /?costo^/ and the computer pre- 
dicts for boresight. This then results in a conical scan. The 
square-law detector or an AGC voltage feeds the ADC 
which is then multiplied by in- and oul-of-phase scan 
sinusoids and is integrated (summed) for one scan. A gain 
h is applied and a teletype printout occurs. If the loop is 
closed, an angular position correction to boresight is made 
and tracking commences. When the system is in closed 
loop operation, the teletype prints the accumulated total 
of corrections. The phase shift 0 is necessary due to phase 



TO SDS 910 
COMPUTER 


Fig. 18. ConicaUtcan hardwara 
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Fig. 19. Conlc«l-tcan toftwara 


lag in the physical antenna and the AGC loop or detector 
filter. The value of 0 is determined experimentally so that 
the two channels decouple and work independently. 

It was found that the primary source of phase lag is 
actually the antenna itself. A circular scan requires 
sinusoidal acceleration in each axis, and the antenna servo 
cannot achieve zero steady state position error. For a 28-s 
scan, B needs to be about - 30 deg and for a 58-s scan has 
to be about - 15 deg in order to decouple the two axes. 
The value of the scan radius R is chosen so that the 


Tabl« 2. Primary claaaaa of um 


Purpose 

Scan 

radius 

Crossover 
loss, dB 

Response 

Maximize gain 

Small 

<0A 

Slow 

Look at dynamics 
(e.g., repeatable 
errors and wind 
effects) 

Large 

2-3 

Fast 


crossover loss is acceptable for the task at hand. For the 
experimental program, each scan had a 2-s halt to allow 
all calculations to be completed. It was not convenient to 
remove this halt, but little error was so introduced because 
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Fig. 20. Croaaovar Iota vs scan radius 
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Tabto 3. Tracking error tatta 


Teit 

Source 

F requency 

Scan radius, deg 

Crossover loss, dB 

Scan period, s 

No. of scans 

rms spread, deg 

Comments 

A 

3C123 

S-band 

0.010 

0.06 

28 

8 

0.00015 

Limited data 

B 

3CI23 

S-band 

0.010 

0.06 

28 

30 

0.00035 


C 

3CI23 

S-band 

0.020 

0.2S 

18 

10 

0.00015 

Antenna jerky, probably 
too short a scan period. 
Limited data. 

D 

3C123 

X-band 

0.002 

0.033 

28 

5 

0.0002 

Little data; however, 
manual hunting could 
only pick up 0.02 dB. 

E 

Pioneer 9 

S-band 

0.040 

1.0 

28 

6 

0.00032 

Limited data. 

F 

Pioneer 9 

S-band 

0.040 

1.0 

28 

8 

0.00035 

Limited data. 

G 

Pioneer 9 

S-band 

0.013 

0.1 

28 

18 

0.0004 


H 

Pioneer 9 

S-band 

0.013 

0.1 

28 

17 

0.00035 



the antenna tended to coast through the halt. The value of 
G is “large” for fast response and “small” for long averag- 
ing times. Two primary classes of use exist as shown in 
Table 2. For ease of relating scan radius to crossover loss, 
Fig. 20 is convenient for S-band (?:i^2300 MHz) and X- 
band («s8400 MHz) frequencies. 


C. Experimental Results 

Two types of data will now be shown — pull-in 
transients and rms error statistics. In Fig. 21 is shown the 
pull-in transient on the radio source 3C123 with the 
S-band beam when the initial error is approximately 
— 0.020 deg in declination. This source gives about 22 
kelvins on the 64-m-diameter antenna at S-band. In Fig. 
22 is shown a pull-in transient on the Pioneer 9 spacecraft 
when initial errors were sizable in both coordinates and 
the initial condition represented a gain loss of 4 dB. 

In Table 3 are shown the results of a series of tests 
performed on October 25, 1973, using HA/DEC scan- 
ning. The crossover loss is taken from Fig. 20. The 
number of scans shown is the number which was observed 
after any initial transient subsided. The rms spread is 



Fig. 21. Conical-scan puli-in on radio source 3C123 at 8-band 


taken as one-sixth of the peak-to-peak variation observed. 
In most cases HA and DEC spreads were comparable. If 
they were not, they were averaged. These data are quite 
crude and obviously are not compensated for any system- 
atic errors. However, the rms spreads arc remarkably 
consistent and may indicate a lower limit to the attainable 
accuracy. In Table 3, only tests B, G and H have a large 
enough number of scans so the rms spread is of value. To 
get more information on the validity of this rms spread, 
refer to Table 4. Tests A, B and C in Table 3 are all on 
3C123. For all times during A, B and C when the system 
was tracking in steady state operation, the overall peak-to- 
peak variation of absolute tracking error vs predicts was 
tabulated. One-sixth of this range is estimated as the rms 
error. The same procedure was applied to the block of 
tests E, F, G and H, Considering that the peak-to-peak 
range observed certainly includes systematic errors of the 
antenna due to the long duration, the rms variation in 
Table 4 probably is a very conservative upper bound on 
the random component of tracking error. It is of interest 
to note that the resolution of the encoders on the master 
equatorial is about 0.0003 deg and their accuracy is on the 
order of 0.001 to 0,002 deg. These data indicate a prob- 
able rms tracking accuracy for moderate strength signals 
like 3C123 and Pioneer 9 of better than a thousandth of a 
degree. This will be discussed further in Section VI-E. 


Table 4. Overall rme variation 


Tests 

Time span 

Observed rms Variation, deg 
Hour angle Declination 

Comments 

A, B,C 

1 h, 13 min 

0.00096 

0.00052 

Begun 1 h after 
meridian 





transit. 

E, F, G, H 

1 h, 59 min 

0.00073 

0.00056 

Straddled 

meridian 

transit. 
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D. Operational Manual Mode 

As the Mariner Venus/Mercury 1973 (Mariner 10) en^ 
counters approached, it was decided to make operational 
for X-band tracking a portion of the experimental conical- 
scan system. The antenna was put into a 58-s scan with a 
radius of 0.(X)4 deg having a crossover loss ^0A4 dB. The 
operator could then average the open loop printouts to 
decide when a correction should be manually entered. 
Alternatively, he could watch the strip chart recorder 
shown in Fig. 18. As boresighl would drift off the 
spacecraft, a sinusoidal pattern would appear on the strip 
chart and the operator could search with HA and DEC 
corrections to null it out. This mode of operation has 
continued through to 1976. 

E. Tracking Data From Mariner Venus/ Mercury 
1973 

During the spring 1974 tracking period of Mariner 10 
considerable tracking time was accumulated in the opera- 
tional manual mode of conical-scan tracking, as described 
above. Each time the operator zeroed the errors, an entry 
was made in a log. In order to get an idea of the potential 
of conical scan (when eventually made fully automatic) 


X 

o 


-10 


-20 h 


-30 


-40 
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▲ 


A A ^ m m 
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• *A 
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■ A 
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TIME (GMT), h 
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A FEB. 13, 1974 


1800 


1900 


Fig. 23. Mariner 10 offseU, February 11-13, 1974 


these log sheets have been reduced to give corrections in 
AZ/EL as well as HA/DEC coordinates and various 
plots were made to try to uncover correlations. 


First, considering repeatability, in Fig. 23 is shown the 
set of raw manual corrections vs time (GMT) for 3 con- 
secutive days. The agreement is within about 0.003 deg for 
the worst case. In Figs. 24 and 25 are shown further 3-day 
plots. Consecutive-day data would be expected to show 
repeatable daily variations, as they do. The worst case 
discrepancy is 0.(X)5 deg for a 1 sigma statistic of less than 
0.(X)1 deg. This is impressive. The consistency indicates 
that a fully automatic conical-scan system will almost 
certainly achieve the rms error of less than 0.001 seen in 
Section VI-C. 


In Fig, 26, elevation offsets vs azimuth are shown for 
March 23-25, 1974. Consistency is again shown, as well as 
a distinct hint that the ionospheric refraction correction is 



Fig. 24. Mariner 10 hour angle offeete, March 23-25, 1974 
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Fig. 27. Marlfiar 10 aiiimith offaala va alavation, Fabruary 4, 1974 


in error for low elevation angles, i.e., when AZ is less than 
140 deg and greater than 240 deg. (A new refraction 
correction program is being prepared at this time.) 

Many plots were made of AAZ vs AZ, AAZ vs EL, AEL 
vs AZ and AEL vs AZ. The curves obtained were rea- 
sonably consistent for consecutive days and/or weeks. 
However, over the roughly four months’ period for which 
much data was taken, only the AAZ vs EL plots gave a 
consistent pattern over the whole time. There are 24 days 
in this period for which a “complete” track was made — 
e.g., low elevation to high elevation to low elevation. 
Every plot of AAZ vs EL showed a characters .tic sideways 
“U”. Figs. 27-31 are typical of this group. The arrow 
indicates time through the track. No reasonable explana- 
tion is offered for this characteristic shape; it may be a 
twisting hysteresis arising in the rising/falling antenna 
structure, but the continuous conical scanning would 
seem to break up such a hysteresis effect. The sides of the 
U became farther apart as summer approached. This is 
probably related to the secant correction for higher eleva- 
tion angles. The U shape may be refraction-connected. In 



Fig. 26. Mariner 10 alavation oftsata va azimuth, March 23^2$, 
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any event it shows the extent to which conical scanning 
can correct for anomalies, whatever their source. 

The data from the logs were subjected to a statistical 
correlation procedure to see if wind effects could correlate 
with any observed offsets. The results were generally 
inconclusive, although on windy days the data seemed to 
scatter more. With a fully automatic conical-scan system, 
the response time can be faster so wind effects can be 
more reliably observed. The problem of correlating sys- 
tematic errors with this antenna is made difficult by the 
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Fig. 29. Mariner 10 azimuth oftaata va alavation, Mareh 9, 1974 
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V. 









ELEVATION, deg 

FI 9 . 30. Mariner 10 ailmulti offaals vt alavatloii, April 1. 1974 



ELEVATION, deg 

FI 9 . 31. Marinar 10 azimuth otfaata va alavatlon, April 11, 1974 

fact that the antenna itself is an AZ/EL mount, while the 
Master Equatorial is of HA/ DEC type. 

F. Comparison of Experimental and Theoretical 
Accuracy 

We are able to make comparisons of the experimental 
and theoretical accuracies for Tests B, G and H of Table 
3 and for the Mariner 10 data. The parameters of these 
tests are given in Table 5. The value of t for B, G and H 
was taken from observation during pull-in, while for 
Mariner 10 it represents a reasonable guess as to the 
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Tablo 5. Tost parsmcisrt 


Parameter 

TestB 

TestG,H 

Manner 10 

R 

0.010 deg 

0.013 deg 

0.004 deg 

W 

0.140 deg 

0.140 deg 

0.038 deg 

T 

75 s 

200 s 

300 s 


20 K 

20 K 

20 K 

L 

22 K 


— 

B 

10 Mm 

— 

— 

Ps 

— 

- 159 dBmW 

- 144 dBmW 


Tabit 6 . Exparimant va thaory 


Test 

Theoretical rms, 
10- ’deg 

Observed rms, 
lO-’deg 

Overall rms, 
10“^ deg 

B 

7.8xlO‘<^ 

0.35 

0.74 

G.H 

1.3 

0.35, 0.4 

0.65 

Mariner 10 

0.045 

<1 

<1 


operator's visual averaging time. In Table 6 the theoretical 
and experimental rms errors are compared. The theoreti- 
cal values ignore gain variations by letting be zero 

in Eqs. (160) and (176). The rms error for Mariner 10 is 
deduced in Section VI-E. The overall rms column is an 
average of the HA and DEC values in Table 4. 

Some interesting conclusions may be drawn from these 
data. First, the observed rms errors are all about the same 
for three vastly differing conditions with greatly different 
theoretical accuracies. This indicates that encoder errors 
and/or gain fluctuations are the primary limitations on 
tracking accuracy. Second, some of the error observed 
might be systematic in nature and the rms error observed 
may not be entirely random. Third, the theoretical value 
of rms error for radio source tracking will be very small 
due primarily to the \/B term in Eq. (160). This small 
value will not be achieved with the present system. 

VII. Conclusions 

A thorough analysis of a conical-scan system suitable 
for large antennas has been made. Details of an experi- 
mental implementation were given and tracking data were 
presented. Tracking accuracies of 0.001 deg rms seem to 
be readily obtainable. 
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Appendix 

Calculation of the Statistics of and 

The means of and must be zero because n(t) was assumed to be white. 
From Eq. (20), 


<Var[Afi'>]' 


n p 

I 


- f^f\('l-'2)S>n“m'lS‘nw«'2‘^'l‘/'2 (A-1) 

•'0 •'0 

where E[ ] denotes expectation and R„{t) is the autocorrelation function of n{t). 
From Eq. (13), 

R.ih- h) - [ Top + Tsg{ R ) ]'«5( - 1,) ( A-2) 

since a Fourier transform relation exists between S„{f) and R„(r). Putting Eq. (A-2) 
in Eg. (A-1) and integrating, we have 

ol-C^[Top-^T,g{R)YBP/2 (A-3) 

Similar manipulation shows that Var[A^‘’^] is also 
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